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ABSTRACT 

We report an investigation of cosmological parameters based on the measurements of anisotropy in the 
cosmic microwave background radiation (CMB) made by Acbar. We use the Acbar data in concert 
with other recent CMB measurements to derive Bayesian estimates of parameters in inflation-motivated 
adiabatic cold dark matter models. We apply a series of additional cosmological constraints on the 
shape and amplitude of the density power spectrum, the Hubble parameter and from supernovae to 
further refine our parameter estimates. Previous estimates of parameters are confirmed, with sensitive 
measurements of the power spectrum now ranging from ^ ~ 3 to 2800. Comparing individual best model 
fits, we find that the addition of Ha as a parameter dramatically improves the fits. We also use the 
high-^ data of Acbar, along with similar data from CBI and BIMA, to investigate potential secondary 
anisotropics from the Sunyaev-Zcldovich effect. We show that the results from the three experiments are 
consistent under this interpretation, and use the data, combined and individually, to estimate as from 
the Sunyaev-Zeldovich component. 

Subject headings: cosmic microwave background — cosmology: parameters 



1. INTRODUCTION 

Anisotropics in the cosmic microwave background radi- 
ation (CMB) are caused by density and temperature fluc- 
tuations in the early universe, when radiation decoupled 
from matter {z ~ 1100). For models in which these initial 
fluctuations are of a Gaussian random nature, the informa- 
tion carried by the CMB anisotropics is completely char- 
acterized by their angular power spectrum as a function 
of Legendre polynomial index £. The comparison between 
measured CMB angular power spectra and theoretical pre- 
dictions can be used to rule out entire classes of cosmolog- 
ical models as well as to estimate the values of a variety of 
cosmological parameters within a given family of models. 

Acoustic oscillations that occurred before the universe 
became largely neutral (z > 1100) produced a harmonic 
series of broad peaks in the power spectrum. For flat mod- 
els, the first peak lies near £ ~ 220. Superposed upon the 
subsequent peak/dip structure is an overall diminishment 
at higher £ because of photon diffusion and the finite thick- 
ness of the last scattering surface. Above £ ^ 1500, the 
entire power spectrum is expected to be strongly atten- 
uated by these effects; this decay at high £ is called the 
damping tail. 

CMB anisotropics have been measured over a wide range 
of £ with high significance. Large angular scale observa- 
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tions with the COBE-DMR instrument produced the first 
convincing detection of CMB anisotropy (Bennett et al. 
1992). At intermediate angular scales, a variety of groups 
have made high signal-to-noise measurements of the first 
peak, near £ ~ 200 (de Bernardis et al. 2000; Hanany et al. 
2000; Halverson et al. 2001; Scott et al. 2002). There is 
also very good evidence for additional harmonic features at 
higher £ (Halverson et al. 2001; Ruhl et al. 2002). Finally, 
at high £ the expected damping tail has been found (Pear- 
son et al. 2002), along with some evidence for greater 
power aX £ > 2000 than is expected for primary CMB 
anisotropics (Mason et al. 2002). 

These measurements provide strong support for the 
inflation- motivated family of cold dark matter models with 
adiabatic initial density perturbations. Improved mea- 
surements of the power spectrum will put these models 
to a more stringent test. The Acbar results reported in 
Kuo et al. (2002, hereafter Paper I) provide the strongest 
constraint to date on the angular power spectrum from 
1000 < £ < 2800, in the damping tail region. In this paper 
we combine these new Acbar results with other published 
CMB results in order to investigate constraints on cosmo- 
logical parameters in the adiabatic A-CDM model space. 

The Acbar results provide detections of power in the 
high-£ region where the effects of secondary anisotropics 
become relevant. In this paper we also combine our high-^ 
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data with previous results from the CBI and BIMA inter- 
ferometers operating at 30 GHz in an attempt to quan- 
tify the contribution from the Sunyaev-Zcldovich Effect 
(SZE) to the power spectrum. Observations such as these, 
at multiple frequencies and over a range of angular scales, 
are essential to the separation of contributions from the 
primary anisotropies and the SZE. 

2. THE INSTRUMENT 

A brief description of the Acbar instrument is given 
in Paper I, with a more complete treatment in Runyan 
et al. (2002). We give only the most relevant details here. 
Acbar is a 16 pixel bolometer array installed on the 2.1m 
Viper telescope at the South Pole, with an angular resolu- 
tion of ~ 5' at 150 GHz. We report here on results derived 
from a subset of those detectors, with bands centered at 
150 GHz. In the first season of observations (2001) the 
array had four such 150 GHz detectors, while in the sec- 
ond season (2002) there were eight. Improvements to the 
receiver and telescope between seasons led to greater sensi- 
tivity and improved pointing reconstruction in the second 
season. 

The Acbar power spectrum used in this work is de- 
rived from observations of two fields covering a total of 
24 deg^ of sky The fields are selected to contain a bright 
guiding quasar near the center, which is removed prior to 
power spectrum estimation. The absolute calibration of 
the instrument is derived from observations of Venus and 
Mars, and has an uncertainty of 10%. The beam profiles 
are derived from images of the guiding quasars and thus 
include any drifts or uncertainties in pointing reconstruc- 
tion. The systematic uncertainty in the beam width, also 
needed here for use in parameter estimation, is 3%. 

3. THE POWER SPECTRUM 

The results of Paper I provide 14 "band powers" with 
effective centers between £ = 190 and 2500. Figure 1 shows 
these results, along with a selection of other current mea- 
surements for context and comparison. The points plotted 
in Figure 1 lie at the maxima Cs of the joint likelihood 
distribution C{{Cb) for each band power (Cb), with the 
vortical bars showing errors derived from the curvature of 
the likelihood at that maximum, 

dCndCB' ■ ^' 

The errors shown are Gaussian ones, [J^~^^^]bb- For 
Acbar, the curvature matrix determined from the orig- 
inal bands has been rotated into the diagonal frame, with 

— 1/2 

eigenvalues /b. Figure 1 shows fg for the errors. For 
the other experiments shown, the results are not rotated 
to the diagonal frame. 

Of course the likelihood curves for the individual band- 
powers are not Gaussians. An offset lognormal distribu- 
tion (Bond et al. 2000) has been shown to be an accurate 
representation of the likelihood curves, and is used in all of 
our parameter estimations. This distribution is Gaussian 
in the variables 

ZB=lli{CB+XB) . (2) 

Here, Xb represents a noise contribution to the band 
power. For AcBAR, xb was treated as a parameter de- 
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termined by fitting the lognormal distribution to the like- 
lihood, as described in in Paper I. Cb, /b^^^i and xb are 
given in Table 2 of that paper. 

4. COMPARISON WITH THEORY 

For the standard cosmological model with adiabatic ini- 
tial density perturbations, the CMB angular power spec- 
trum can be readily calculated as a function of input cos- 
mological parameters y. These theoretical predictions for 
Ce can be made individually for each i, while our measure- 
ment is over bands of finite width, characterized by the 
window functions (fse- Given the theoretically predicted 
Ci, the predicted band power in one of our bands is 

Cb = l{ipBtCi)/I{ipBi) , (3) 
where Ce = i{i + 1)^/277 and 

The window functions ipBe give the response of the 
bands to power at each £. Numerical tabulations of the 
window functions are available on the Acbar experiment 
public website^^. They have been rotated into the diago- 
nal frame. 

The likelihood that the data would result from the cos- 
mology described by y is given by 

ln£(C) =ln£(C)- 

IY.^Zb- Zb)T^^UZb' - Zb') , 

BB' 

where Zb = ln{CB + xb) is the value of the lognormal 
parameter at the position of maximum likelihood Cb and 
J^BB' the curvature matrix transformed into the lognor- 
mal variables, 

^bI = (Cb + xb):F^bI' (c'b + x'b) • (6) 
Provided with the maximum likelihood band powers Cb, 

(C) 

the lognormal offsets xb , the curvature matrix ^, , and 
the window functions ipBt, the likeHhood of the parameter 
set y, given the data, can be computed. 

Our set y consists of seven cosmological parameters: 
rifc = (1 - ^tot),^tv,LOcdni,LOb,ns,Tc,\TiCw. The total en- 
ergy density of the universe in units of critical density 
Q,tot is linked to the global curvature of space: negatively 
curved for Q.tot < 1, positively curved for ^tot > 1 and 
flat for fltot = 1- The total energy density has 3 con- 
stituents: vacuum (Oa), matter and rclativistic particles. 
The rclativistic energy density is currently negligible. The 
matter density is split into two types, baryonic matter 
(fib = ujb/h'^), which interacts with electromagnetic radia- 
tion, and cold dark matter (flcdm = i^'cdm/h^), which does 
not. The total matter density is denoted Qrn = ^b + ^cdm- 
The amplitude of the CMB power spectrum at £ = 10, 
InCio gives the overall amplitude of the primordial fluc- 
tuations. This quantity is well-constrained by the COBE- 
DMR observations (Bennett et al. 1996). The full COBE- 
DMR power spectrum as described in Bond et al. (2000) 
is included in all of our parameter analyses. The spectral 
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index of primordial density perturbations, ng, parameter- 
izes the variation in the fluctuation power as a function of 
length scale; = 1 corresponds to scale invariance. 

The universe reionized at some point between decou- 
pling and the present. After reionization, CMB photons 
scatter further, tc is the Compton optical depth (from 
decoupling to present) due to such scattering. High tc di- 
minishes CMB power by a factor of exp[— 2tc] over most 
of the £ range, though not in the DMR range. 

Many more parameters than our basic 7 may be needed 
to completely describe inflationary models. These include 
the gravity- wave induced tensor amplitude and tilt, vari- 
ations of tilt with wavenumber, relativistic particle den- 
sities, more complex dynamics associated with the dark 
energy Ct\, etc. For example, a gravity-wave induced com- 
ponent is expected in the largest class of inflation models, 
and has impact on the spectrum in the £ < 100 region. Al- 
though far from the region where AcBAR has its impact, it 
can affect the amplitude, tilt, and Compton depth. If one 
constrains the tilts and amplitudes of the tensor compo- 
nent to those motivated by simple inflation models, there 
is little change in the other parameters. 

It is possible that secondary anisotropies, such as the 
Sunyaev-Zeldovich effect investigated later in this paper, 
could contribute significant power to the highest £ band of 
our measurement. However, the detection of anisotropy in 
that band is only l.lcr, or 0.9a above the best-fit model 
primary CMB angular power spectrum. Thus, for the pur- 
pose of cosmological parameter estimation from the pri- 
mary CMB signal, we can safely ignore the effects of po- 
tential SZE contamination. 

To derive estimates of cosmological parameters, we com- 
pare our data with the primary CMB power spectra pre- 
dicted by combinations of cosmological parameters y. In 
this comparison, we vary InCio continuously, while the 
other parameters take the discrete values listed in Table 1. 
An angular power spectrum is generated for each set of 
discrete parameter values, forming a grid. Given a CMB 
dataset (e.g., Acbar or a combination of various measured 
power spectra), the likelihood is then calculated for each 
point on the seven dimensional grid. 

To compute the likelihood that a particular parameter 
X has a value xq, the seven dimensional grid of likeli- 
hoods is integrated over the other six parameters, holding 
the parameter of interest fixed at xq- This method, known 
as marginalization, involves calculating 

L{X = xo) = J 5{X - xo)Pprioriy)C{y)dy , (7) 

where d{x) is the usual delta function and the prior 
Pprior{y) is discusscd bclow. 

For each model on the parameter grid, along with the 
overall amplitude parameter Cio, we continuously vary 
the beam widths ahi and calibrations Ai for each ex- 
periment i about their estimated values a^i and unity, 
to take into account the uncertainties in the respective 
measurements. We approximate the beam-uncertainty 
and calibration-uncertainty "prior" probabilities by Gaus- 
sians in A((T(,i)^ and InAf, respectively. The modifi- 
cation to the bandpower as a result of the uncertainty 
A(CTbi)^ is modeled by exp[-((^ -|- ^f)B^{(Jbif], with 

((^+5)^)5 =I{^Be{£+^ f}/I{^Be)- The ovcraU impact 
of the calibration and beam uncertainties is that the com- 



bination \nCiQ + J2i lii^i - ((^+ 5)^)BA(cr6i)^ is adjusted 
for each grid parameter combination to give the best fit 

with errors. Marginalization over the continuous parame- 
ters is done by calculating a Fisher matrix, and assuming a 
Gaussian distribution in the posterior distributions in the 
{InCio, In^i , ^(iTfci)^} variables. Marginalization over the 
grid parameters is done by discrete integration. 

In addition to the parameters given in Table 1. we ex- 
amine four "derived" parameters, which are functions of 
our basic 7 variables, {to, h, al, and Feff}. Here h is the 
value of the Hubble expansion parameter, Hq = 100 h 
km/s/Mpc and is given by 

The age of the universe is 

, 9.778 2x^ 
io(Gyr) « — - / dx , (9) 

where we have dropped the minor effects of relativistic par- 
ticles here, but not when we do our actual comparisons. 
The variance in the (linear) density fluctuation spectrum 
on the scale of clusters of galaxies {8h~^ Mpc) is The 
shape of the linear density power spectrum is described by 
the parameter 

Feff = F {n, - l)/2, where 

r = flrnh e-(^''(l+f^m v^)-0.06) _ ^^^^^ 

Note that FeS ~ ^Imh with corrections due to baryon 
density and spectral tilt rig over the region probed by 
large scale structure observations that approximate the 
dominant dependences. The derived parameters are not 
marginalized over, and do not define dimensions in the 
parameter grid because they are determined given a set 
of parameters y. The derived parameters are calculated 
at each point on the model grid. We can make Bayesian 
estimates of those parameters using equation (7) with X 
being the derived parameter of interest. 

We can also use the derived parameters to account for 
other cosmological information. Estimates of these param- 
eters from other, non-CMB, observations can be included 
as "prior constraints" . Equation (7) is cast in a form that 
gives a likelihood as a function of our parameters, and 

priors Pprioriv)- 

All our analyses include the loose prior implicitly im- 
posed by the edges of the database given in Table 1, and 
two other very weak priors that are generally accepted by 
most cosmologists: models are restricted to those where 
the current age of the universe is to > lOGyr, and where 
the matter density flm is greater than 0.1. 

We add to these constraints a series of additional priors: 

• weak-h : 0.45 < h < 0.90. This is a tophat restric- 
tion, designed to allow the CMB data, rather than 
arguable priors, to drive the results. 

• LSS : We employ two constraints based on large 
scale structure (LSS) observations, on the combina- 
tions (TgfiJ^'''^ and Fcfi'. For both priors, we adopt a 
Gaussian distribution convolved with a top-hat dis- 
tribution, characterized by the parameters agft'^^ = 

0.47+;o2;l:;i8 and Fcff = 0.21+:o3;l-°s, where the first 
error gives the 1-a point of the Gaussian distribution. 
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and the second error gives the extent of the tophat 
distribution. Our basic philosophy is to adopt priors 
that are not overly restrictive since the LSS data is 
still improving. The motivation for the choice and 
the discussion of the LSS data is given in Bond et al. 
(2002a). The Feff distribution encompasses recent 
results from the 2dF and SDSS surveys. The as 
distribution encompasses results from recent weak 
lensing surveys. It also covers many of the cluster 
abundance determinations using X-ray temperature 
and other cluster data. 

• LSS(Zo?i)— (Tg): There are currently a few cluster 
abundance estimations that point to values of as 
that are lower than the weak lensing estimates. Al- 
though our standard LSS prior takes most of these 
variations into account by its spread, we have also 
tested the effect of shifting the entire ag-distribution 
downward by 15%, to agfi^^^ = 0-^Ot-°Ql't;ll. We 
keep the F prior the same. 

• HST-h : We strengthen our h prior, based on 
the Hubble Space Telescope Key Project measure- 
ment (Freedman et al. 2001) of the Hubble constant, 
h = 0.72 ± 0.08. This is a Gaussian prior with the 
stated error as the 1-a points. 

• Strong data : We combine the LSS prior given 
above with the HST-/i prior, and add a constraint 
(in the Oa vs. ^Itot plane) based on surveys of the 
brightness vs. redshift relation of Type la Super- 
novae (Perlmutter & Riess 1999; Riess et al. 1998). 

• Flat : Inflation models generally predict a flat 
geometry, and recent evidence supports this (eg 
de Bernardis et al. 2000; Halverson et al. 2001). For 
the converted, we investigate the effects of holding 
Otot equal to one. 

5. CONSTRAINTS ON COSMOLOGICAL PARAMETERS 
FROM CMB SPECTRA 

We apply these methods to three combinations of CMB 
data, using a series of priors for each case. By doing so, 
we can investigate the power of adding Acbar data to 
the current cosmological mix, as well as the dependence of 
results on the strength and nature of the applied priors. 

The full COBE-DMR power spectrum of Bond et al. 
(1998) is included in all our analyses. The first combina- 
tion of CMB data, Acbar+DMR, investigates the poten- 
tial for using the DMR low-^ anchor with the damping-tail 
measurement of AcBAR as an independent check on pre- 
vious CMB-based cosmological parameter estimation. 

Figure 2 shows marginalized likelihood curves for 
AcBAR and DMR only, using the wcak-h and weak-ft,+flat 
priors. These are displayed for each of two sets of AcBAR 
data, the first being the full spectrum (bands 1-14), the 
second being bands 2-14. The curves show the parameter 
estimates are very prior-dependent, and also depend upon 
whether we include the lowest i band data, despite the 
large errors in that band. These variations are a result 
of AcBAR not pinning down the peak/dip structure well. 
Parameter estimates arc then relying more on the specific 
shape of the damping tail. 



The physics of damping is well known and clearly de- 
pends upon the cosmological parameter characterizing the 
strength of the viscous and diffusive couplings, LUh. It 
might then be thought that by using only Acbar and 
DMR we could get a strong constraint on this parameter. 
However, when we take account of all of the influences that 
determine the damping scale in ^-space, it is found to be 
relatively insensitive (see e.g., Sicvers ct al. (2002) for a dis- 
cussion). Once some information is given on the peak/dip 
structure that helps to pin the parameters, AcBAR im- 
proves the quality of the determinations by virtue of its 
small error bars in the damping tail region. We therefore 
proceed by including CMB data in the £ range between 
DMR and Acbar, over the first three acoustic peaks. 

Many CMB experiments have made sensitive measure- 
ments of the power spectrum in the region of those peaks; 
we choose here to combine the low-^ measurement of DMR 
with a set of recent higher-i? observations, comprised of 
Archeops (Benoit et al. 2002), Boomerang (Ruhl et al. 
2002), DASI (Halverson et al. 2001), MAXIMA (Hanany 
et al. 2000), and VSA (Scott et al. 2002). We also include 
the recent high-^ results of CBI (Pearson et al. 2002). We 
give the label "Others" to this aggregate set of data and 
first investigate the parameter extraction that can be done 
with these measurements, sans Acbar. Finally, we inves- 
tigate the improvements made when adding AcBAR to the 
mix, in a combination we label "ACBAR+Others". 

Two plots of one-dimensional marginalized likelihood 
curves, for "Others" and "AcBAR-|- Others" , are given in 
Figures 3 and 4 respectively. Here we see that results are 
generally stable to the application of priors - that is, the 
application of a prior may narrow the result, but does not 
move it outside the range imposed by other priors. 

It is remarkable how many parameters are well con- 
strained in Figures 3 and 4. Using only a weak-h prior, four 
of the five parameters shown {ilk, ^bh^, ^cdmh^, and Ug) 
are well localized. The detection of becomes stronger 
as stronger priors are applied. 

We have investigated the effects of dropping the first 
Acbar bin on the results shown in Figure 4, and unlike 
the AcBAR-fDMR results, the effects of this action are 
negligible. This is not surprising, given that the first peak 
information in "Others" and "Others-|-AcBAR" is driven 
by the other measurements. 

The likelihood curves shown in these figures can be used 
to find confidence intervals on each parameter for each 
prior. We can proceed in a similar way to constrain the 
values of other parameters, such as h and the current age 
of the universe (to). Table 2 gives the median values and 
1-a limits (16% and 84% integrals of the likelihood) for 
a set of cosmological parameters. For tc, the 95% confi- 
dence upper limit is given. A comparison of the values in 
the Table, or of Figures 3 and 4, shows that the impact 
of adding AcBAR to this mix is not dramatic. There are 
modest improvements in some parameter estimates, most 
notably Oa and Qcdmh^- 

Not immediately obvious from the figures is the im- 
proved rejection of fl^ = models. One measure of this 
is the improvement in the 3-a lower limit on f^A, found 
by integrating the likelihood curve. This 3-a limit im- 
proves from Qa > 0.086 for "Others" to Qa > 0.136 
for "AcBAR-|-Others" ; that is, the probability of having 
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a lower value of than these is 0.14%. This depends 
upon the specific range in h we are allowing in our weak 
prior. 

The source of this rejection of models near ~ can 
be illustrated by examining the of the aggregate dataset 
to the best-fit models, in both the = and the "free" 
cases. We find, for the "AcBAR+Others" dataset (con- 
sisting of 116 band powers) ~ 140 and x^ = 160 for the 
best-fit "free JIa" and f^A = models, respectively. Thus, 
while both models plotted in Figure 1 appear reasonable 
to the eye, the fit is significantly better for the f^A = 0.5 
model. 

Unfortunately, calculating the effective number of de- 
grees of freedom in this x^ is not straightforward. Taking 
beam and calibration uncertainties for each observation as 
a total of 16 parameters in the fit added to the 7 cosmolog- 
ical parameters, we know the effective degrees of freedom 
lies in the range 93 < dof < 116. Adopting dof=100 as a 
reasonable estimate, the probability of finding x^ > 140 
and x^ > 160 are P> = 0.0051 and 0.00013 respectively. 
We caution the reader against strict interpretation of these 
statistically high x'^ in the face of this very heterogeneous 
data set. Instead, we note the significant improvement in 
enabled by the addition of a single parameter in the fit. 

Interestingly, the x^ difference of the best-fit models is 
roughly consistent with the ratio of the marginalized like- 
lihoods (for the weak-ft, prior) at the peak of the likelihood 
curve near f^A = 0.7 and its level at Qa = 0. Exact corre- 
spondence would be expected if our parameter likelihood 
function had a Gaussian form. 

One would think that the information ACBAR adds at 
high I would significantly improve the determinations of 
Wfo because of viscous damping and rig just because of the 
increased £-baseline. This is clearly not the case. We have 
discussed the damping tail issue already. The near degen- 
eracies due to correlations among certain parameter com- 
binations imply that increased data does not necessarily 
lead to increased precision on the cosmological parameters. 
It is likely that the lack of improvement in these projec- 
tions to individual parameters is due to degeneracies in 
the parameter space. 

Such degeneracies are well known and have been dis- 
cussed at length in the literature; see Efstathiou & Bond 
(1999) for an extensive treatment of some of the most per- 
nicious of these. For example, one of these (the "geometric 
degeneracy" ), leads to nearly identical angular power spec- 
tra for particular combinations of (I^a, ^^fc, f^m, /i), while 
leaving and Wcdm fixed. The breaking of this geometric 
degeneracy is the reason estimates of JIa and, to a lesser 
extent ^cdm^ ■, improve so dramatically as stronger priors 
on h are applied. Other less exact degeneracies, such as 
one between amplitude, and rc, lead to similar broad- 
enings of these projected likelihood curves. 

By choosing our database parameters well (e.g., using 
the physical densities w;, and u^cdm rather than the densi- 
ties relative to critical) we have minimized the effects of 
some potential degeneracies. By exploring the parameter 
eigenmodes we can escape the limitations of the canonical 
parameters and determine the true power of any combi- 
nation of datasets. This process is in fact quite familiar 
to most cosmologists; the Type la supernovae results are 
often considered as limiting the parameter "eigenmode" 



r^A-^^M- 

Table 3 lists the best-determined five (of seven total) 

eigenmodes for the "Others" and "AcBAR-|-Others" analy- 
ses. The coefficients describing those modes, and the errors 
on the eigenvalues, are determined by ensemble averages 
of the likelihood derivatives over the parameter database. 
We introduce the probability-weighted ensemble average 
of a parameter j/^, 

(ya) = y yaPprior{y)C{y)dy , (11) 

and the probability weighted ensemble average of the dif- 
ferentials 

{5va 5yb) = j 5ya5ybPprior{y)C{ij)dy , (12) 

where 5ya = ya - {ya}- 

The eigenmodes themselves (^fc) and their errors (cr^) 
are given by the application of a rotation matrix R to the 
parameter vector y, 

Cfe = ^ RakSya , (13) 

a 

{SyaSyb) = Rak{4fRbk ■ (14) 

k 

In this eigenmode analysis, we have used the fractional 
deviations Sub/ (tOb) and Sucdm/ {^cdm) as parameters, 
rather than Ub and Ucdm, to set their deviation magnitudes 
on more equal footing with those of the other parameters. 

Inspection of the table shows that while the eigenmodes 
for "Others" and "AcBAR-f- Others" are not identical, they 
are very similar. In most cases they are dominated by con- 
tributions from one or two cosmological parameters, but 
in all cases there are significant components from several 
parameters. The table shows much more clearly the im- 
pact of adding Acbar to the dataset; all of the eigenvalue 
uncertainties improve, but the greatest improvement is in 
the fifth eigenmode, which becomes the fourth one when 
Acbar is added. It is dominated by the cosmological con- 
stant. 

6. CONSTRAINTS ON CTg FROM THE SUNYAEV-ZELDOVICH 

EFFECT 

The Acbar results provide the first data above £ = 2000 
at 150 GHz. The recent CBI Deep field results (Mason 
et al. 2002), at 30 GHz, have indicated a possible excess 
over the expected primary anisotropy signal at ^ > 2000. 
The most promising candidate for the source of the excess 
is the Sunyaev-Zeldovich Effect (SZE) due to the scatter- 
ing of CMB off hot electrons in the intra-cluster medium 
(see Birkinshaw (1999) for a recent review). The CBI 
results have been interpreted in the context of the SZE 
with tentative constraints being obtained on the value of 
as (Bond et al. 2002a; Komatsu & Seljak 2002; Holder 
2002). The BIMA array (Dawson et al. 2002), operating 
at 30 GHz, has also reported detection of power at higher 
i which also has been attributed to the SZE by Komatsu 
& Seljak (2002). 

Parameter fitting using sccoiidary effects such as the 
SZE must be approached with caution. Both numerical 
and analytical predictions for the SZE power spectrum 
suffer from a number of uncertainties. The results of dif- 
ferent simulations, although in general agreement, show 
significant differences in both the amplitude and shape of 
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the predicted spectrum. Analytical models suffer from un- 
certainties inherent in modeling the profile of the clusters. 
In addition, cooling and heating effects in the clusters are 
not yet well understood and most simulations and analyt- 
ical models do not take these effects into account. Simu- 
lations have shown that for modest deviations about the 
concordance ACDM model, the SZE angular power spec- 
trum scales as ^ (f7{,ft,)^(7|. Despite uncertainties in 
the physics of cluster models, especially the role of energy 
injection, and the relatively large errors on the observa- 
tions, this very strong dependence of the SZE spectrum 
on as enables the derivation of significant constraints on 
as from current data. 

We choose to model the SZE using two angular power 
spectrum templates. The first is obtained from large 
Smoothed Particle Hydrodynamics (SPH) simulations of 
ACDM (Bond et al. 2002b; Wadsley et al. 2002). The 
second is obtained from an analytical model (sec Zhang 
et al. (2002); Bond et al. (2002a) for details). Both tem- 
plates were scaled to a fiducial value of Ofe/i = 0.035. 

Although the power in the primary spectrum is falling 
rapidly compared to the rising contribution of the SZE at 
£ > 2000, interpretation of the low noise Acbar band 
powers around the cross-over region is sensitive to the 
contribution of the primary signal together with the sec- 
ondary. Rather than consider a full range of parameter 
space we select a simple model for the primary spectrum, 
a best fit flat model for the "ACBAR+Others" data com- 
bination (Ob = 0.047, ^cdm = 0.253, riA = 0.7, h = 0.69, 
ris = 0.975 and tc = 0.2). The primary model was nor- 
malized with the best fit amplitude obtained from the 
fits. However, uncertainties in the model parameters affect 
the overall amplitude of the primary spectrum at scales 
£ pa 2000 where the primary and secondary signals are 
comparable. We chose to parameterize the freedom in the 
primary and secondary spectra by two effective parame- 
ters, an amplitude in the primary power spectrum 
and a scaling factor for the SZE af^. Uncertainty in the 
primary parameter 172^ represents the uncertainties in a 
number of dominant effects given by the combination of 
parameters: as, tc, and shown in section 5. Most 

significantly, it also incorporates the effect of systematic 
uncertainties such as an overall calibration and beam un- 
certainty in the data. The secondary amplitude parameter 
(jf ^ describes the scaling of the SZE spectrum and can be 
related to as via af^ « {nbh/0m5)°-'^^as. 

We select points with £ > 1500 for fitting and use the 
offset lognormal approximation as described in Section 3. 
The target model is now Zg = \n{CB + g^C^ + xb) where 
the primary and SZE spectra have been filtered by the ap- 
propriate window functions for each band power and Qi, 
is the frequency dependent scaling of the SZE (a factor 
of ~ 4 in power lower at 150GHz compared to 30GHz). 
Top-hat windows are used to model the BIMA band pow- 
ers. The primary model is scaled by gr^^ over the range 
(0.1, 1.8) and the secondary model is scaled by af^ over 
the range (0.5, 1.4). At each point in the grid, we calculate 
the quantity —2 ln(L) for the model. 

To assign a realistic uncertainty to the amplitude of the 
primary models, we add a Gaussian prior in the amplitude 
with a width of 20% RMS. This is chosen to reflect 
the uncertainty in the amplitude of the best fit models 



obtained in the parameter fits described in section 5. As 
an example, the Acbar data fixes the amplitude of our 
template model with an RMS of 17% while the CBI Deep 
field data fixes the same model with an RMS of 21%. We 
check that our results are robust to a change in the width 
of the prior by fitting with a 10% and 40% RMS width. 
The marginalized, best fit value, and upper errors for crg^ 
reported below changes by only 0.1% while the lower error 
changes by 1% on average. 

In Fig. 5 we show contour plots of the grids in the 
(glir, ff^) plane. We subtract the value at the mini- 
mum from the grids and the (2.3, 6.17, 11.8) contours give 
an indication of where the 1, 2, 3 a levels would fall if the 
likelihoods were Gaussian. We show results for both tem- 
plate models used in the analysis. The contours show that 
the data is only weakly dependent on the amplitude 
with a slight movement to higher uf ^ for low as ex- 
pected for Acbar and CBI data. This reflects the fact that 
the only bands sensitive to the primary amplitude are the 
lowest £ bands included from the Acbar and CBI data. 
The high £ BIMA data is fully degenerate in the amplitude 
of the primary, as expected, and the contours in the di- 
rection simply reflect the Gaussian prior. It does however 
give a strong lower bound although this depends strongly 
on the lognormality of the BIMA band power distribution. 
The xb offsets arc not available for BIMA and are set to 
xb = although we also show results where no lognormal 
transformation was applied to the BIMA points (Table 4). 
Acbar and CBI provide strong upper bounds. A com- 
bination of the three data sets give strong constraints in 
o"!^ with the contours showing a slight tilt in the expected 
direction with respect to ^3^. 

One major obstacle in using a secondary effect such as 
the SZE to fit parameters is the non-Gaussian nature of the 
signal from non-linear structures such as clusters. In gen- 
eral the non-Gaussianity will increase the sample variance 
of the underlying signal. Treating the data as Gaussian 
can therefore result in an overestimate of the significance 
of the constraint. This effect has been investigated us- 
ing numerical simulations of the SZE (White et al. 2002; 
Zhang et al. 2002) and also by calculating the contribu- 
tion to the covariance by the fourth order, trispcctrum 
term T«/ (Cooray 2001; Komatsu & Seljak 2002). In gen- 
eral the sample variance is found to be a factor ^ 3 higher 
than the Gaussian equivalent with some dependence on £ 
and also on the width of the bands being considered {£-£' 
correlations are also altered by the non-Gaussianity). In 
order to include this effect in the errors and correlations 
we scale the inverse Fisher matrix of the band powers as 

f— 1 = fngTT- 1 fUg Cir^ 
•'BB' — J B •'BB'Jb' 1 \^^> 

where is the scaling factor. Here we chose such 
that the sample variance component of the error is a factor 
of 3 larger than the expected Gaussian case which can be 
approximately evaluated as 



V/skyAB(24 + l) ■ 

Fig. 6 shows the effect of the correction on {q2%:,(^t'^)- In 
general the correction changes from experiment to exper- 
iment and from band to band due to the different sample 
variance component in the uncertainty of each band. It is 
therefore important to include this effect in the analysis 
since it can have a substantial effect on the structure of 
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the contours as opposed to a simple rescaling of the con- 
fidence Umits. We find that the correction has significant 
effects on the allowed region, particularly at the 2 to 3-a 
level. 

To obtain best estimates on the value of crf^, we 
marginalize over the g^i: direction to recover the one di- 
mensional likelihood in crf^. We show the resulting like- 
lihoods in Fig. 7 for both analytical and SPH templates. 
Both results include the non-Gaussian correction discussed 
above. The 2-a region and median values shown as error- 
bars are obtained by calculating the 2.5%, 97.5% and 50% 
integrals of the likelihoods respectively. The results for 
both templates with 2-(t error estimates arc summarized 
in Table 4. We find that fitting with the SPH template 
results in values for crf^ about 6% higher than when the 
analytical model is used. This effect is due to the SPH 
model having a lower amplitude at larger scales than the 
analytical model (Fig. 8). 

In Fig. 8 we show the primary and SZE templates used, 
scaled to the best fit value of crf^ = 0.98 (at 30GHz and 
150GHz) for the analytical case shown in Fig. 7. The total 
primary-t-SZE (analytical) model is also shown together 
with the SPH template scaled to the same parameters. 
We see that the method can obtain a good fit between the 
data and primary-t-SZE model at both observing frequen- 
cies. Fig. 8 also shows the non-Gaussian corrections to 
each band power error. 

It is important to note that if there is a non-negligible 
SZE component to the observed power, it may affect the 
parameter fits which assume only a primary contribution. 
However, the relative contribution to the ACBAR band 
powers is very small, only approximately 15% in the last 
three bands. We do not expect this to have any signif- 
icant impact on the parameter estimates derived in this 
work. As future observations increase the accuracy in this 
region of the spectrum, a fully consistent approach to pa- 
rameter fitting will have to be adopted. Such an approach 
would simultaneously account for primary and secondary 
anisotropy, fitting for all parameters. 

7. CONCLUSIONS 

The AcBAR data, the most sensitive to date in the 
damping-tail region, are in good agreement with predic- 
tions of flat A-CDM models with adiabatic initial pertur- 
bations. Considering the AcBAR data together with other 
recent CMB results, we find that the addition of a single 
parameter to the model, CI a, dramatically improves the 
best fit, with an improvement of Ax^ = 20 upon adding 
that one parameter. Using very weak cosmological priors 
(Om > 0.1, age > lOGyr, 0.45 <h< 0.90), the 3-a lower 
limit on the cosmological constant rises to > 0.136 
upon including the Acbar data. 

We find that the addition of Acbar data to the cur- 
rent CMB set does not lead to substantial improvements 
in the 1-a estimates of the canonical cosmological param- 



eters. However, in an eigenmode analysis, the addition of 
Acbar data does improve the rotated parameter uncer- 
tainties, indicating that in this case the lack of improved 
errors on the pure cosmological parameters is probably 
dominated by degeneracies between those parameters. 

We fit a SZE component to the data using Acbar and 
other measurements at high-^. Our estimates for the value 
of the effective quantity (t|^ ~ {flbh /0. 035)^ "^^ as show an 
improvement over previous estimates using only CBI and 
BIMA observations. Although CBI and Acbar alone do 
not provide lower bounds on trf^, the combination of the 
two observations results in a detection greater than 3-a 
which is independent of the BIMA lower bound. 

Our fits show that at a fiducial value of ilbh = 0.035 the 
central values for as are consistently higher than other es- 
timates obtained using cluster data, weak lensing surveys 
and primary CMB observations (see Bond et al. (2002a) 
for a recent survey of estimates). However, the results 
overlap with most other estimates at the 2-a level. As 
an indication, our LSS prior, a smoothed tophat on the 
combination as^rn"^, translates roughly into a smoothed 
tophat of CTg — 0.921q']^5 at flm = 0.3. Any statistical 
inconsistency therefore, appears to be mild. Furthermore, 
systematic uncertainties in the estimates have not yet been 
taken into account. The difference displayed by the numer- 
ical and analytical based results of about 10% is indicative 
of the agreement between the two methods for predicting 
the SZE power spectrum (Komatsu & Seljak 2002). In 
addition, entropy injection may have a significant effect 
on the SZE power spectrum. These effects would change 
the shape of the SZE template and impact directly on our 
determination of as- Nevertheless, we conclude that the 
in the context of the phenomenological models adopted in 
this work, the data is consistent with a SZE component at 
as values near the high end of independent estimates. 
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Fig. 1. — Top Panel: The Acbar CMB power spectrum, Ce = £{£ + l)Ci/{2n), plotted over a vacuum energy dominated 
[flk = -0.05, = 0.5, ujcdm = 0.12, uJb = 0.02, Hq = 50, rc = 0.025, Us = 0.925, amplitude VCio = l-H x IO'^^Tcmb] model 
(black thin line) and a CDM dominated [Qk ~ 0.05, Qa ~ 0, ujcdm ~ 0.22, cut = 0.02, Ho — 50, tc — 0, = 0.925, amplitude 
VCio = 1.34 X IO'^Tcmb] model (green thick line). These are the best-fit models, for A and A-free models respectively, found 
during the Acbar-|- Others parameter estimation described in the text, with the weak-h prior. Bottom Panel: The top panel with 
the addition of power spectra from several other experiments. Both models appear to be reasonable fits to the data, with the 
= 0.5 model statistically being the better of the two. 
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Table 1 

Parameter Values Used in Grid 
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Note. — Grid point values for the six cosmological parameters that are varied discretely. The grid is not 
spaced evenly for Wcdm, wj, na, and tc; points are purposefully more concentrated in regions in which the 
likelihood (found from previous datasets) is high. We only calculate models on this grid which have > 0.1; 
this, along with the edges of each parameter range, forms an implicit prior in our analysis. 
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Parameter Estimates and Errors 
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Note. — Parameter estimates and errors for several prior combinations with and without Acbar,. Errors are quoted at 1 cr (16% 
and 84% points of the integral of the likelihood), except for tc where the 95% upper- limit is given. The various priors are described 
in the text. The top block lists results found with and without the inclusion of Acbar data, which shows the small improvements 
found upon adding Acbar to the mix. The bottom block shows the effect of applying stronger priors on the AcBAR+Others dataset, 
which naturally leads to dramatic improvements on the parameter estimates. The difference between the LSS and LSS(low-cr8) 
priors (discussed in the text) does lead to several slight shifts, smaller than the 1-cr errors. 
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Table 3 

ElGENMODES 



Eigenmode 


Error 


TO 










Us 




"Others" — Archcops^- 


Boomcrang+CBI+DASI 


+DMR+Maxima-|-VSA 






1 


0.012 


-0.233 


0.082 


-0.084 


0.111 


-0.724 


0.604 


0.174 


2 


0.017 


-0.210 


-0.223 


-0.007 


-0.238 


0.607 


0.676 


0.155 


3 


0.046 


0.432 


0.236 


-0.751 


0.146 


0.111 


0.231 


-0.327 


4 


0.091 


0.107 


0.700 


0.521 


0.315 


0.196 


0.244 


-0.171 


5 


0.126 


0.421 


0.007 


0.260 


-0.750 


-0.231 


0.157 


-0.340 


"ACBAR+Others" 
















1 


0.010 


-0.140 


0.133 


-0.051 


0.194 


-0.889 


0.344 


0.123 


2 


0.015 


-0.264 


-0.155 


-0.026 


-0.238 


0.320 


0.843 


0.189 


3 


0.043 


0.489 


0.339 


-0.673 


0.070 


0.081 


0.259 


-0.338 


4 


0.063 


0.368 


0.279 


0.362 


-0.780 


-0.198 


0.050 


-0.078 


5 


0.088 


0.012 


0.630 


0.541 


0.451 


0.223 


0.206 


-0.123 



Note. — Eigenmodes and the uncertainties on their determination, for the 
"Others" analysis (top set) and the "AcBAR-|-Others" analysis (bottom set). 
Only the top five (of seven) are listed. The first column labels the eigenmodes 
in rank order of uncertainty; these uncertainties, and the eigenvectors chosen, 

arc derived from weak-ft prior ensemble averages over the database as described 
in the text. The second column lists the uncertainty on each eigenmode, while 
columns 3-9 list the coefficients of the eigenmode rotation matrix Rak, applied 
to the basis set of parameters labelled at the top of the columns. There is 
significant improvement, especially for the fifth eigenmode, upon adding Acbar 
to the dataset; note that it becomes the fourth eigenmode of "AcBAR-fOthers" , 
and shows a factor of two improvement. 



Table 4 
Estimates 



Data SPH Analytic 



ACBAR 


1.04°; 


17 

48 


1.12°: 


.18 

.55 


ACBAR+CBI+BIMA 


1.04^; 


11 
21 


0.98" 


12 
21 


ACBAR-l-CBI-l-BIMAt 


i.o4«: 


12 
42 


0.96o 


.15 
.41 



Note. — The results for the phenomeno- 
logical fits to a SZE component. The table 
shows crf^ estimates with statistical 2-a er- 
rors, for both SPH and analytical models. 
Non-Gaussian corrections axe included, 
t No lognormal transformation applied to 
the BIMA band power. 
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Fig. 2. — Likelihood curves for Acbar+DMR with the weak and uieafc+flat priors. Each prior case is plotted with and without 
the first ACBAR band (centered on I = 187) included in the analysis. In the upper-right panel, the 1 to 3-(t contours are shown for 
the 2D r^fc-r^A likelihood with the weak prior and bands 1-14 (blue) and bands 2-14 (red). The thick black lines define Q,m = 
and Q,m = 1 and the dotted black line defines Q.m ~ 0.5. The yellow contours are the 1, 2 and 3-a levels of constraints based on 
Type la Supernovae. The lack of stability of the curves (for u>i, in particular) indicates that Acbar-|-DMR alone is not sufficient 
for robust parameter estimation. 
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Fig. 3.— Likelihood curves for Archeops+B98+CBI+DASI+DMR+MAXIMA+VSA ("Others") with the weak, lueafc+LSS, 
HST-^, and strong data priors. The ^k-^h contours are shown for the weak (blue) and strong data (red) cases. The yellow 
contours are the 1, 2 and 3-(t levels of constraints based on Type la Supernovae. CMB estimates of Q.k, cocdm, and uit are stable 
with sensible behavior as additional priors are employed. 
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Fig. 4.— Likelihood curves for Acbar +Arclieops+B98+CBH-DASI+DMR+MAXIMA+VSA (Acbar+ "Others") with the 
weak, weak+LSS, HST-/i, and strong data priors. The fifc-ilA contours are shown for the weak (blue) and strong data (red) cases. 
The yellow contours are the 1, 2 and 3-a levels of constraints based on Type la Supernovae. The positions and widths of these 
curves do not differ significantly from those in Figure 3 despite the addition of the low noise Acbar data through the damping 
tail. A comparison of the curves here containing the LSS prior {weak+LSS and strong data) with those derived using the a lower 
estimate (discussed in the text) shows only small changes. The most noticeable changes are an upward shift in the lower tail on 
SIa, and a broader and higher- value peak. Table 2 gives numerical estimates of these parameters, derived by integration of 
these curves. 
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^sz 

Fig. 5. — 1, 2, and 3-a contours for various combinations of data sets in the (g2A')'^l^) plane. The left panels show fits 
obtained using the SPH template and the right panels are the equivalent for the analytical model. Both CBI and Acbar data 
points constrain the upper values of crf^ but not the lower values as they are sensitive to the amplitude of the primary spectrum. 
The higher-^ BIMA observations are insensitive to the primary component and therefore provide a strong lower bound. The 
combination of the three datasets (bottom row - dotted, black contours) show a strong detection of the SZ component. The 
dashed parallel lines show the width of the Gaussian prior imposed on q2%- We use a lognormal distribution for the BIMA band 
powers. 



16 



Goldstein et al. 



model = SPH sim 

1 1 1 1 1 1 1 1 1 1 1 1 1 r 




Fig. 6. — The plot shows the effect of adding a correction for the increased sample variance of the data due to the non-Gaussian 
scatter. Both contours are for the combination Acbar+CBI+BIMA. The effect is significant, particularly at the 3-a level. 
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Fig. 7. — Marginalized likelihoods for the Acbar+CBI+BIMA combination. The results of fits using both the SPH (red, dashed) 
and analytical (blue, solid) models are shown. Both cases include non-Gaussian corrections. The 2-cr upper and lower bounds and 
median values were obtained by computing the 2.5%, 97.5%, and 50% integrals of the likelihoods respectively. We find that the 
SPH model prefers slightly higher values for crf^. This is due to the fact that the SPH spectra show less power than the analytical 
models on large scales. 
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1000 104 

Fig. 8. — Best fit primary+SZ model amplitudes. The top panel shows the SZ model scaled to an observing ficciucncy of 30GHz 
which corresponds to the BIMA and CBI observations (first four rightmost points and two leftmost points respectively). The lower 
panel is for 150GHz, corresponding to the Acbar results. The solid lines show the analytical model at each frequency and the 
total primaxy+SZ. The dotted line shows the fiducial spectrum used to model the primary contribution at ^ > 1000. The dashed 
line shows the SPH spectrum used as a template shape. The SZ contributions are scaled to a value of uf^ = 0.98 to show the 
amplitude at the best fit value obtained from the ACBAR+CBI+BIMA combination using the analytic model as template (with 
non-Gaussian corrections). The extended errorbars for each point show the corrections due to the increase in sample variance 
expected from the non-Gaussian signal. 



